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OPERATOR SCALING STABLE RANDOM FIELDS 



HERMINE BIERME, MARK M. MEERSCHAERT, AND HANS-PETER SCHEFFLER 

Abstract. A scalar valued random field {X{x)}^^^d is called operator-scaling if 
for some dx d matrix E with positive real parts of the eigenvalues and some H > 
we have 

{X{c^x)}^jzud ^= {c^X(x)}^gKd for all c> 0, 

where denotes equality of all finite-dimensional marginal distributions. We 
present a moving average and a harmonizable representation of stable operator 
scaling random fields by utilizing so called i?-homogeneous functions (p, satisfying 
(p{c^x) = c(p{x). These fields also have stationary increments and are stochasti- 
cally continuous. In the Gaussian case critical Holder-exponents and the Hausdorff- 
dimension of the sample paths are also obtained. 



1. Introduction 

A scalar valued random field {X{x)}^^^d is called operator-scaling if for some dx d 
matrix E with positive real parts of the eigenvalues and some H > we have 

(1.1) {X(c^x) W {c^X(x)K.gM^ foralloO, 

where denotes equality of all finite-dimensional marginal distributions. As usual 
= exp(_E'logc) where exp(A) = J2T=o^ matrix exponential. Note that if 

E = I, the identity matrix, then (11.11) is just the well-known self-similarity property 
{X (ex)} rcfz-^d ^= {c^X(x)}^g]Rd where one usually calls H the Hurst-index. See [TT] 
for an overview of self-similar processes in the one-dimensional case d = 1. Self- 
similar processes are used in various fields of applications such as internet traffic 
modelling [2H], ground water modelling and mathematical finance, just to mention a 
few. Various examples can be found for instance in the books [TH] and [T]. A very 
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important class of such fields or processes are the fractional stable fields and especially 
the Levy fractional Brownian field. 

These fields have different definitions which are usually not equivalent. More pre- 
cisely, for < a < 2 let Za{dy) be an independently scattered symmetric a-stable 
(SaS) random measure on with Lebesgue control measure A'' (see [21] P- 121). 
For < H < 1 one defines the moving average representation by 

(1.2) Xh{x)= [ (\\x-y\\^~'^/^-\\y\\^-''/'^)z^{dy). 

For Waid^ a complex isotropic SaS random measure with Lebesgue control measure 
the harmonizable representation is given by 

(1.3) Xnix) = Re / (e^<-'«> - \)\\i\\-''-''l'^WMi)- 

See for a comprehensive introduction to random integrals with respect to stable 
measures. It follows from basic properties that {XH{cx)}x£Rd ^= {c^X(x)}^gKd as 
well as {XH{cx)}x^^d ^= {c^ XH{x)}x^^d. Moreover, both processes have stationary 
increments, that is for any h we have {XH{x + h)—XH{h)}x(z^d ^= {XH{x)}x(zRd 
and similarly for {XH{x)}x^^d. Furthermore both fields are isotropic, that is 
{XH{Ax)}x(z;^d ^= {X{x)}x(z^d for any orthogonal matrix A. It is worth mention- 
ing that if a < 2 the fields {XH{x)}xeRd and {XHix)}^^'' defined in ()1.2|1 and ()1.3|) . 
respectively, are usually different. See [21], Theorem 7.7.4 for the one-dimensional 
case. However, in the Gaussian case a = 2, by computing the covariance function of 
the fields, it follows that {XH{x)}x^T^d and {XH{x)}^^^d have the same law up to a 
multiplicative constant and known as the Levy fractional Brownian field. 

Certain applications (see, e.g., [3 |H1 ED] and references therein) require that the 
random field is anisotropic and satisfies a scaling relation. This scaling relation should 
have different Hurst indices in different directions and these directions should not 
necessarily be orthogonal. In the Gaussian case a prominent example of an anisotropic 
random field is the fractional Brownian sheet {BH{x)}xeM.d defined as follows: Let 
< ifj- < 1 for j = 1, . . . , (i and set H = {Hi, . . . , Ha). Define 

Bh{x) = [ f\[\x,- u,\^^-'/' - \u,f^-'/'] Z,{du). 
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See HH EH and the literature cited there for more information on these fields. 
Then, if we set E = dia,g{H^^ , . . . , H^^), it follows by a simple computation that 
{BH{c^x)}x(z]^d ^= {c-B//(x)}^giKd so {BH{x)}x(z^d is operator scaling in the sense of 
(jl.ip . However, {i?//(x)}^g]Rd does not have stationary increments. 

The purpose of this paper is to define two different classes of operator scaling stable 
random fields (OSSRF) and analyze their basic properties. We present a moving 
average representation as well as a harmonizable representation. Our constructions 
are based on a class of E'-homogeneous functions ip -.W^ ^ [0, oo) where if is positive 
on M'^\{0} and ip{c^x) = c(p{x) for all x G M'^ and c > 0. Such functions were studied 
in detail in [12], Chapter 5. It will turn out that the harmonizable representation 
allows more fiexibility in the class of possible functions ip in contrast to the moving 
average representation which is more restrictive. However, in both cases the OSSRFs 
satisfy (ILljl . have stationary increments and are continuous in probability. In the 
Gaussian case a = 2 we show that there exists modifications of these fields which are 
almost surely Holder-continuous of certain indices and we compute the Hausdorff- 
dimension of the graph. 

This paper is organized as follows: In section 2 we introduce the class of E- 
homogeneous functions, derive some basic properties and provide important examples. 
In section 3 we define and analyze a moving average representation of OSSRFs. Sec- 
tion 4 is devoted to the harmonizable representation and its properties. Finally, in 
the Gaussian case a = 2, we analyze the sample path properties of both the moving 
average and the harmonizable representation of OSSRFs. 

2. E-HOMOGENEOUS FUNCTIONS 

Let E he a real d x d matrix with positive real parts of the eigenvalues < ai < 
a2 ■ ■ ■ < ttp for p < d. Let us denote F = M*^ \ {0}. It follows from Lemma 6.1.5 of 
[TH] that there exists a norm || ■ ||o on M'^ such that for the unit sphere 5*0 = {x G M'^ : 
||x||o = 1} the mapping \i/ : (0, oo) x So — F, \i/(r, 6^) = r^9 is a homeomorphism. 
Moreover for any x G F the function t i-^ ||^^a^||o is strictly increasing. Hence we 
can write any x G F uniquely as x = r(x)^/(x) for some radial part r(x) > and 
some direction /(x) G Sq such that x i— r(x) and x i— /(x) are continuous. Observe 
that 5*0 = {x G M'^ : r(x) = 1} is compact. Moreover we know that r(x) ^ oo as 
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X ^ oo and t{x) ^ as x ^ 0. Hence we can extend r(-) continuously by setting 
r(0) = 0. Note that further t{—x) = t{x) and l{—x) = —l{x). The following result 
gives bounds on the growth rate of t{x) in terms of the real parts of the eigenvalues 
of 

Lemma 2.1. For any (small) S > there exist constants Ci, . . . , C4 > such that 
for all \\x\\o < 1 or all t{x) < 1, 

and, for all \\x\\o > 1 or all r(x) > 1, 

C3||x||;/"^~'<r(x) <C4||x||;/''^+'. 

Proof. We will only prove the first two inequalities. It follows from Theorem 2.2.4 of 
[ini that for any 5' > we have t'^^"''' ||t-^6'||o as t ^ 00 uniformly in ||6'||o = 1. 
Hence \\t~^\\o := sup^g^^ ll^~'^6'||o < Ct~"i+'^' for alH > 1 and some constant C > 0. 
Equivalently ||s-^||o < Cs"^^"^' for all s < 1. Since ||x||o = ||r(x)-^/(x) |||^ < ||t(x)-^||o < 
Cr(x)"^~'^' we get t{x) > Ci||x||q^"^''''^, for 6 = ^^^r:^ ~ ^7 if ll^^llo < 1 which is 
equivalent to r(x) < 1. 

Similarly we know that, for any 5' > 0, t~'^p~^' \\t^6\\o ^ as t — 00 uniformly in 
ll^llo = 1. Therefore ||t^||o < 0^""+^' for alH > 1 or equivalently ||s-^||o < Cs^""^^' 
for all s < 1. But x = t{x)^1{x) and l{x) = t{x)"^x. Thus, 1 < ||''"(3;)~^||Q||a^||Q and 
ll^^llo > C~W{xfp+^' for all ||x||o < 1. Hence r(x) < CaHxH J^"^"^ for 5 = ^ - ^ 
and ||x||o < 1. The proof is complete. □ 

The following results generalize some of the results in [12], Chapter l.A to our more 
general case of exponents E. 

Lemma 2.2. There exists a constant K >1 such that for all x, ?/ G M'^ we have 

rix + y) < K{t{x) +r(y)). 

Proof. Observe that the set G = {(x, G M'^ x M"^ : r(x) + T{y) = 1} is bounded 
by Lemma l2T] and closed by continuity of r. Hence G is a compact set. Thus the 
continuous function (x, t-^ r(x + y) assumes a finite maximum K on G. Since 
'S'o X {0} C G, we have K > 1. Given any x, ?/ G M'^ both not equal to zero we set 
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s = (r(x) + r(y))~\ Then, with t{c^x) = ct{x) it follows that 

r(x + y) = s^V(s-^(x + y)^ = s~^t{{s^x) + (s^y)). 

But (s^x, s^y) E G since r (s-^x) + r (s-^y) = s (r(x) + 'r(i/)) = 1. Therefore, 

t{x + y)< Ks~^ = K[t{x) + T{y)) 

and the proof is complete. □ 

Now let q = trace(i?) and observe that by multivariable change of variables we 
have X'^{c^{B)) = c''A'^(i?) for all Borel sets B C M^, c > 0, which can be written as 
d{c^x) = d^dx. Let B{r,x) = {y eM.^ : T{y — x) < r} denote the ball of radius r > 
around x G M'^. Then it is easy to see that B{r, x) = x + B{r, 0) = x + B{1, 0) and 
hence A'^(i?(r, x)) = r'^A'^(-B(l, 0)). The following proposition provides an integration 
in polar coordinates formula. 

Proposition 2.3. There exists a unique finite Radon measure a on Sq such that for 
all f G L^(M'^, dx) we have 

poo p 

f{x)dx= / / f{r^9)a{d9)r''~^dr. 
Jo Jsq 

The proof of Proposition 12.31 is based on the following. 

Lemma 2.4. If f : T ^ C is continuous and f{r^x) = r~'^f{x) for all r > and 
X G r, then there exists a constant fif such that for all g E ((0, oo), r~^dr^ we have 

f°° dr 

f{x)g{r{x))dx = fif g{r)—. 



r 



Proof. Let Lj : (0, oo) ^ C be defined as 

{fi^ / ^^ fix) dx if r > 1 

-Il<r{x)<r-^f(^)dx ifr<l. 

Since / is continuous on F, from dominated convergence, Lf is continuous on (0, 1) U 
(1,+cx)). But \'^{B{r,0)) = r«A'^(5(l,0)) implies that A'^({x G R'^ : r(x) = r}) = 0, 
and it follows that Lf is also continuous at point 1 and thus on (0, +oo). Moreover, 
for any r > we have Lf{r~^) = —Lf{r). When rs > 1 with r, s > a change of 
variables yields 

Lf{rs)= I f{x)dx= I fis''y)sUy= [ f{y)dy. 

Jl<T{x)<rs J l<T{s^y)<rs Js ^<T{y)<.r 
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Let US assume for instance that 1 < s^^ < r. Then, by continuity of Lj, 

/ f{y) dy= [ f{y) dy- [ f{y) dy = Lf{r) - Lf{s-'). 

It follows using Lf{s^^) = — that 

(2.1) Lfirs) = Lfir) + Lfis). 

Similarly we show that (12.111 holds for s^^ < I < r and < r < 1 and thus for all 
rs > 1. Using again the fact that Lf{r~^) = —Lf{r)^ for all r > 0, H^'I.l^ is valid for all 
r,s > 0. By continuity of Lj it follows that Lf{r) = Lf{e) logr. We set /i/ = Lf{e). 
If g{r) = l]a,b](r) for some < a < 6 we get 



f{x)g{T{x))dx = / f{x)dx = Lf{b) - Lf{a) 

'a<T(x)<b 



f°° dr 

fXf (log b - log a) = fXf g{r)—. 

Jo 



The general result follows by taking linear combinations and limits of these functions 
in the standard way. □ 

Proof of Proposition \2.f\ When / G C{Sq) define / on F by /(x) = r(x)~'^/(/(x)) . 
The function / satisfies the hypothesis of Lemma EH If / > then /i/ = Lj{e) = 
Ii<r{x)<e ^i^yfiK^)) dx > 0. Moreover = afij, = fij+fig and the mapping 
/ (— s> /ij is continuous. Hence this mapping is a positive linear functional on C{So). 
Therefore there exists a Radon measure a on 5*0 such that /i/ = Jg^ f{0) a{d6). 

gi G Cc((0, oo)) we get from applying Lemma EH with / and g{r) = r'^gi{r) that 



fil{x))gi{T{x))dx = / f{x)T{xygi{T{x))dx 
= fif I giiry'^dr 



/ f{e)a{de)g^{ry-Ur. 
J So 

Since linear combinations of functions of the form f{l{x))gi{r{x)) are dense in 
L^(M'^, dx) the result follows. □ 

Corollary 2.5. Let /3 G M and suppose f : ^ C is measurable such that \f{x)\ = 
0(r(x)^). If P > —q then f is integrable nearO, and if P < —q then f is integrable 
near infinity. 
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We are now in position to define the class of -E-homogeneous functions and a 
important subclass needed in the moving average representation of OSSRFs. Let E 
he a d X d matrix as above such that < ai < ■ ■ ■ < and for x G F let {t{x), l{x)) 
be the polar coordinates associated to that is a; = t{x)^1{x). 

Definition 2.6. Let (p : ^ C he any function. We say that (p is E-homogeneous 
if (p{c^x) = ap{x) for all c > and x eT. 

It follows that an -E-homogeneous function (f is completely determined by its values 
on So, since (f{x) = i^{t{x)^1{x)) = T{x)ip{l{x)). Observe that ii (p is -E-homogeneous 
and continuous with positive values on F, then 

(2.2) = maxip{9) > and = minv9(^) > 0. 

OgSq 0&Sq 

Moreover by continuity we necessarily have ^{0) = 0. 

Definition 2.7. Let j3 > 0. A continuous function (yj : M°' ^ [0, oo) is called {P,E)- 
admissible, if ip{x) > for all x 7^ and for any < A < B there exists a positive 
constant C > such that, for A < \\y\\ < B, 

r{x) < 1 ^ \^{x + y)- ^{y)\ < Cr{xY . 

Remark 2.8. If a continuous function ip : W^' ^ [0, 00) is positive and Lipschitz on 
F, that is |</?(x) — 'p{y)\ < C\\x — y\\o for x,y eT, then (p is -E')-admissible for all 
/9 < ai by Lemma ITTl 

Remark 2.9. If (f is £;)-admissible then (3 < ai. In fact, if = Vi ® ■ ■ ■ ® Vp is 
the spectral decomposition of M*^ with respect to E (see [11], Chapter 2 for details), 
by restricting the argument of the proof of Lemma 12.11 to the space Vi on can show 
that for any 6 > there exists a constant C > such that t{x) < C\\x\\l^^^ ^ for 
all X E Vi with ||a;||o < 1. Then, if for some fixed nonzero n G Vi we consider the 
function t 1—^ ip{tu) we get for Si = (35 that \(p{tu + su) — 'p{su)\ < for all 

small t and s bounded away from zero and infinity. If one would have /3 > ai, one 
could chose 6 > such that P/ai — 61 > 1 and hence there would exist a constant 
K > such that ip(tu) = K for all t 7^ 0. But since </? is continuous and <^(0) = 
this is impossible. 
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Remark 2.10. In general the exponent E of a homogeneous function ip is not unique. 
It is easy to check that ^{x) — ^ oo as ||x|| oo, and then Theorem 5.2.13 in ^H] 
implies that the set of possible exponents is E + TS{ip) where E is any exponent, 
S{lp) is the set of symmetries of and TS{ip) is the tangent space at the identity. 
Here we say that A is a symmetry of (f if (f{Ax) = (f{x) for all x G M'^. The 
symmetries S{(p) form a Lie group, and the tangent space consists of all derivatives 
x'(0) of smooth curves x{t) on S{ip) for which x(0) = / the identity. For example, if 
if is rotationally invariant then S{(f) is the orthogonal group and TS{(f) is the linear 
space of skew-symmetric matrices. Although exponents are not unique, Theorem 
5.2.14 in [in] shows that every exponent -E of a homogeneous function ip has the 
same real spectrum < ai < ■ ■ • < and induces the same spectral decomposition 
M.'^ = Vi Q) ■ ■ . (BVp, since these structural components describe the growth properties 
of the homogeneous function. In particular, the function r i— > Lp(rx) grows like r^/"^ 
for any nonzero x E Vi, see Section 5.3 in [in] for more details. 

We conclude this section by examples of -E')-admissible, E'-homogeneous func- 
tions ip : M.'^ [0, oo) used in Theorem 13.11 below to define a moving average repre- 
sentation of OSSRFs {X^(x)}^gKd. Let us denote < ., . > the standard inner product 
on and -E* the transpose of any d x (i-matrix E with respect to this inner product. 
The following class of examples is inspired by the log-characteristic function of a full 
operator stable law on M.'^. See [IH] for details. 

Theorem 2.11. Assume E is a real d x d-matrix such that the real parts of the 
eigenvalues satisfy 1/2 < Oi < 02... < ctp for p < d. Assume M{d6) is a finite 
measure on the unit sphere Sq corresponding to E such that 

spanjr-^'^ : r > 0,^ G supp(M)} = R'^. 

Then 

r foo , dr 

if{x)= / {1 -cos{{x, r"" 6))) —M{de) 

J So Jo 

is a continuous, E-homogeneous function such that (p{x) > for all x eT. Moreover 
ip is {P, E) -admissible for (3 < min [ai, if < ^ and P = lifai>l. 

Proof. Let ai > 1/2 denote the smallest real part of the eigenvalues of E. Since 
E and E^ have the same eigenvalues, it follows from Theorem 2.2.4 of ^H] that 
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for any 5 > there exists a constant C > such that ||?"^* 6*110 < Cr"-^~^ for all 
< r < 1 and 9q G Sq. Therefore, from dominated convergence, ip is well-defined 
and continuous on M'^. Moreover we have ^{x) > and ^{x) = implies x = 0. 
A simple change of variable shows that ip{c^x) = ap{x) for all c > and x G Mf^. 
It remains to show that ip is i?)- admissible. Using the trigonometric identity 
cos(a) — cos(6) = —2 sin((a + b)/2) sin((a — b)/2) we have for any x,y eM.'^ that 

{x,r^'e)- 



(2.3) \ipix + y)-ipiy)\<2 



So JO 



'{x + 2y,r^'t 

sml 



sm 



uT 

-^M{d9). 



First, let us assume that ai > 1, then an upper bound of ()2.3|) is given by 

■(x,r^'6')' 



So JG 



sm 



nT 



which is finite because oi > 1, using ||r^*6'||o < Cr°-^~^ for all < r < 1 and 6*0 G 5*0, 
and elementary estimates. Moreover writing x = t{x)^1{x) a change of variables 
yields to 



So Jo 



sm 



^M{de) = 2r{x) [ [ 

^ J So JO 



^a(x),r'''9) 

sm' 



flT 

-M{de), 



which proves that ip is 1-admissible. 

Let us now consider the case where Oi < 1. Choose 6 > small enough. On one 
hand, for r < 1, one can find C > such that 

{x + 2y,r^'e)\ ^,_j{x,r^'e)- 



sm 



sm 



<c(2\\y\ 



+ ||x||o ||x||or 



On the other hand, it follows from Theorem 2.2.4 of |Tn| that one can find C > 
such that ||r^'6'||o < Cr""-^^ for all r > 1 and 9o G 5*0. Thus, for 7 < min (^1,^), 
using I sin(n)| < one can find C > such that 

{x + 2y,r^'9)\ . /{x,r^'9)' 



sm 



sm 



< C\\x\\2r"'''+^^. 



Therefore, by substituting these upper bounds into the right-hand side of (I2.3jl and 
integrating, for some constant C > we have shown that \ip{x + y) — f{y)\ < Cll^^llo 
for all ||x||o < 1 and A < \\y\\o < B. 

Since by Lemma ED ||ic||o < Ct{xY'^~^ for r(x) < 1, the assertion follows with 
/3 = 7(ai-(5). □ 
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The following result gives a constructive description of a large class of continuous, 
admissible i5-homogeneous functions. 

Corollary 2.12. Let 9i,...,9d be any basis ofW^', let < \i < ■ ■ ■ < \d and 
Ci, . . . ,Cd > 0. Choose a d x d matrix E such that E^Oj = XjOj for j = 1, . . . ,d. 
Then for any p > 0, if p < 2Ai the function 

is a continuous E-homogeneous and {f3, E)-admissible function for {3 < min ^Ai, p^^ 
if ^1 < P and [3 = p if Xi> p. 



Proof. First observe that since r^^Oj = r^^Oj it follows that '^{c^x) = C!f{x). Moreover 
ip is continuous. Let B > A > 0, since y i— X]j=i ^j) 1^^'^'' continuous and 

positive on F, by the mean value theorem, for A < \\y\\ < B and < v4/2, one can 
find C > such that 

d d 

(2.4) W{x + y)- ^{y) \ < c| C, | (x + 9,) - C, \ (y, 9,) ^1^^ . 

3=1 i=i 

Hence it remains to show that the right hand side of ()2.4|) is -E)-admissible. Let 
=i7j^0j for suitable 7j > 0, where Se denotes the dirac mass in 9. Let us 
define for a; G M*^, 

^{x) = / (1 - cos((x, r^^/")^^))) ^ M{d9), 

J So Jo 

which is well defined since p < 2Ai. Moreover, by Theorem I2.1H ip is {j3, {1/ p)E)- 
admissible for P < min if Ai < p and /5 = 1 if Ai > p. Let Tp{x) denote 

the radial part with respect to [1/ p)E. Then uniqueness implies that the radial part 
t{x) with respect to E is given by r(x) = Tp^xY^^. Hence if) is -E)-admissible for 
(3 < min ( Ai, p^ ) if Ai < p and /3 = p if Ai > p. 
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Moreover, since r^^^f^^^Oj = r^^^^Oj we get 

m=Y.^, / (l-cos(rVp|(^,^^.)|)) 

= V ^ ( r (1 - cos(s))s-(''/^^)-^ ds) I (x, r/^^ 

This completes the proof. □ 

3. Moving average representation 

In this section we consider a moving average representation of OSSRFs and derive 
its basic properties. We first give sufficient conditions such that the integral repre- 
sentation exists. More precisely, for < a < 2 we consider Za{dy) an independently 
scattered SaS random measure on with Lebesgue control measure A*^. Then we 
define a moving average representation of OSSRFs using the basic fact that a random 
integral j^^ f{y) Za{dy) exists if and only if J^^ \ f {y)\°' dy < oo. 

Throughout this section we fix a real d x d matrix E with < ai < ■ ■ ■ < 
denoting the real parts of the eigenvalues of E. As before, let q = trace(_E'). 

Theorem 3.1. Let /3 > 0. Let (p : M.'^ ^ [0,oo) be an E -homogeneous, {(3,E)- 
admissible function. Then for any < a < 2 and any < H < (3 the random field 

(3.1) X^{x) = [ Ux - y)^-^/" - ^{-yf-'^/A Z^idy) , X G 

exists and is stochastically continuous. 

Proof. Let us recall that X^{x) exists if and only if 



^{X - y)^-"/" - ^{-y)^-l/^ 



dy < oo. 



Let us assume that H G (0,/5). Observe that by (|2.2|) and the fact that ip is E- 
homogeneous, ip{z) < M^t{z) and (p{z) > m^r{z) for all 2; 7^ 0. Fix any x G F. 
Then, 

V9(x - y)^-'^/" - ^(-y)^-'^/"P < c{r{x - yr""-" + r(y)°^-^) . 



12 



HERMINE BIERME, MARK M. MEERSCHAERT, AND HANS-PETER SCHEFFLER 



But for any i? > it follows from Corollary 12.51 that X-(y)<R '^dy < oo if 



H > 0. Moreover, by Lemma IT2l {^ : T{x — y) < R} C {y : T{y) < K{R + t{x))} and 
hence, by a change of variable we obtain using Corollarv 12.51 again that, if H > 

t{x - y)""-'^ dy = 



r(y)<R 



T{yY^~^dy < / T{yY^~^dy < oo. 

r{x-y)<R JT{y)<K(R+T{x)) 



It remains to show that for some R = R{x) > we have 



(3.2) 



T{y)>R 



dy < oo. 



^{x + y)^-"/" - y.(y)^-''/" 

Observe that for r(y) > R, (p{y) > 0, so we can write 

ip(x + y) = ^ {^{y)^ {^{y)~^x + ^{yY^y))) = ip{y)ip {ip{yY^x + ip{yY^y) , 

since ip is -E-homogeneous. Moreover ip {fiyY^v) = 1 since ip is {P,E) admissi- 
ble, one can find C > such that 

|(p ((p(y)-^x + p^iyY^y) - l\ < Cr (^(y)-^x)^ = C^iyY^r {xf . 

Hence by the mean value theorem applied to the function t^"?/" ^gar t = 1, one can 
find Ci > such that 



p}{x + y) 



H~q/a 



v{y) 



H~q/a 



= v^ly)^""/" kiviyY^x + viyYSf-''"' - i| 

< {xf , 

for all T{y) > R, where i? > is chosen sufficiently large so that C(p{yY^T (x)^ < 1/2 
for all r{y) > R. But v?(y)^"'^"^/" < C2r(?/)^-^-9/" and by Corollary EIH we know 
that j^i^y)>p, r{yY^~''^"''^ dy < oo li H < (3. This allows to conclude that r^(x) is 
finite for all x G M"'. Let us now show that is stochastically continuous. Since 
X^p is a SaS field, it follows from Proposition 3.5.1 in [21] that X^p is stochastically 
continuous if and only if, for all Xq G M'', 

<^(xo + x- 1/)^-^/" - <^(xo - y)^-"/" 

By a change a variables, this holds if and only if 

(3.3) r°(x) ^0 as X ^ 0. 

But if is continuous on so 



dy as X ^ 0. 



(^(x - y) 



H-q/a _ 



vi-y) 



H~q/o 



as X ^ 
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for almost every y G W^. Moreover, arguing as above, as soon as r(x) < 1, for suitable 
i? > 0, one can find C > such that 

where 1b (y) denotes the indicator function of a set B. Then ()3.3|) holds using domi- 
nated convergence. 

This concludes the proof. □ 

Corollary 3.2. Under the conditions of Theorem \S.l\ the random field {X^(x)}^.g]gd 
has the following properties: 

(a) operator scaling, that is, for any c > 0, 

(3.4) {X^(c^x)}^gMd =' {c^X^(x)}^.gMd- 

(b) stationary increments, that is, for any h G M'^, 

(3.5) {X^{x + h) - X^{h)}^^^d =' {X^(x)}^gK''- 

Proof. We will only prove part (a). The proof of part (b) is left to the reader. Fix 
any Xi, . . . ,Xm. G M.'^. Then (j3.4^) follows if we can show that for any ti, . . . ,tm G M 
we have 

m m 

^t,X^(c^x,) = c^^t,X^(x,). 
i=i j=i 

By a change of variable together with ip{c^x) = ap{x) and the fact that Za{c^dz) = 
c''/°'Za{dz) we get 

Y,t,X^{c^x,)= / ^tY^(c^a;,-y)^-^/"-v^(-y)^-^/")z„(rfy) 



= c^^t,X^(x,) 

and the proof is complete. □ 
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Remark 3.3. Theorem 13.11 and Corollary 13.21 inchide the following classical isotropic 
random fields as special cases. Assume (f{x) = \\x\\ and E = I, the identity matrix. 
Observe that ip is an ii^-homogeneous, (1, £')-admissible function. Then 

X^i^)= ! {\\x-y\f -'I'' -\\y\f Z^{dy) 

Especially, if a = 2, then {X^{x)}x^-^d is known as the Levy fractional Brownian field. 
Note that in this case, for any < a < 2 equation ()3.4I) reduces to the well-known 
self-similarity property {X^{cx)}x(zKd ^= {c^ X^{x)}x^^d. Moreover our results also 
include the well known one-dimensional case c? = 1 of linear fractional stable motions 
and especially the fractional Brownian motion when a = 2. 



4. Harmonizable representation 

In this section we consider an harmonizable representation of OSSRFs and derive 
its basic properties. We first give necessary and sufficient conditions such that the 
integral representation exists and yields a stochastically continuous field. For < 
a < 2, let Wa{d^) be a complex isotropic SaS random measure with Lebesgue control 
measure (see [21] p. 281). 

Throughout this section we fix a real d x d matrix E with < ai < ■ ■ ■ < Op 
denoting the real parts of the eigenvalues of E. As before, let q = trace(-E). 

Theorem 4.1. Let ip : 'R'^ [0, oo) be a continuous, E^ -homogeneous function such 
that ■ip{x) 7^ for x ^ 0. Then for any < a < 2 the random field 

(4.1) X^x) = Re [ (e*<-'«> - l)V'(0~^"'^" W^id^ , x G M'' 

exists and is stochastically continuous if and only if H E (0,ai). 

Proof. Let us recall that X^{x) exists if and only if 

r°(x):= / |e'<"'«>-l|°V'(0""''"'^^^<+oo. 

Let us assume that H G (0,ai). By integration in polar coordinates for E^ given by 
Proposition 12. 3| 



So 



^^^V(^)""''^''cr(c/^) dr. 
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For 6 E {0,H — ai), by considering the cases r > 1 and < r < 1 separately and 
using the same spectral bounds on the growth of ||r^* || as in the proof of Lemma l2?T| 
one can find C > such that 

- 1 ° < C (1 + min {r'^(''^-^\ 1). 

Moreover, since ifj is continuous with positive values on the sphere 5*0, and hence 
bounded away from zero, 



So 

This allows us to conclude that r^(x) is finite for all x G M*^. Let us show now that 
is stochastically continuous. Since X^, is a SaS field, it is stochastically continuous 
if and only if, for all Xq G M'^, 

((e'<"°+"'«> - 1) - (e'<-°'«> - 1)) ^0 as a; ^ 

that is, equivalently, 



(4.2) 



r^(x) -^0 as X ^ 0. 



It is straightforward to see that ()4.2|) holds for H G (0,ai), using dominated conver- 
gence and the upper bound computed above. 

Conversely, let us assume that X^ exists and that it is stochastically continuous. Let 
us remark that in this case r^(x) exists for all x G M"^ and satisfies, for all A > 

r^(A^x) = A"^r;(x). 

Let us fix any x G M*^, with x 7^ and let us notice that r^(x) 7^ 0. Since X^ is 
stochastically continuous, by ()4.2p 



A"^r;^(x) -^0 as A ^ 0, 

which implies that H > 0. 
Let us now prove that H < ai. 

First case : Assume that ai is an eigenvalue of E. Then there exist 61 G such that 
1 1 6*1 1 1 = 1 and E61 = ai^i. Therefore 



JSo 



aH—q 



with 



\<ei,r^ e>\ = r"^| <ei,e>\ <Cr''\ 
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\ 1/ai 



Then, for r < (§) 

\^i<e„r^'e> _ ^1 _ 2|sin ( "'^^'' ^ ) \ > 2r 



TT 



and hence 



'So 

Since is positive on the sphere Sq, 



\< 9i,9>\''ij{9y"-''(T{d9) > 



So 



and then r^(^i) < +oo imphes that H < ai. 

Second case : Assume that ai is not an eigenvalue of E. Then there exists 6i G M such 
that Ai = «! + ibi and Ai are complex eigenvalues of E. One can find ^1,71 G M*^, 
with ||6'i|| = II71II = 1 such that 

r^9i = r"^ (cos (61 logr)6'i + sin (61 logr) 7i) 

r^7i = r"^ (— sin (61 logr) 6^1 + cos (61 logr)7i) . 

Then it can be shown using the inequality je*"^ — M > 1^1/^ for 1^1 < ^ that a lower 
bound of T^{9i) + r^(7i) is given by 

TT" JQ J So 

Observe that for a, 6 > we have a" + 6° > (a^ + h'^y''^. Therefore 



|2 



a/2 



|<r%,^>|" + |<r%,^>|° > (|< r%,^ >|% |< r%,e >| 

> r"*^! (|<^l,^>|^ + |<7l,0>|2)"^^ 

Then we conclude as in the first case that H < ai. The proof is complete. □ 

Corollary 4.2. Under the conditions of Theorem \4.1[ the random field {X^{x)}x£Rd 
has the following properties: 

(a) operator scaling, that is, for any c > 0, 

(4.3) {X^{c'^x)}^^^d =' {c-^X^(x)}^giRd. 

(b) stationary increments, that is, for any h G M"', 

(4.4) {X^{x + h) - X^(/i)}^gKd ^= {X^(x)}^gKd 
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Proof. Let us recall that by corollary 6.3.2 of for / G L°'{M.'^), the characteristic 
function of the random variable F = Re J^^ f{y)Wa{dy) is given by 

(4.5) E(e^*^) =expf-co|tr / \f{yTdy) where co = ^ T {cos 9f d9. 

Hence, for any Xi,...,Xm G M'^, the finite dimensional characteristic function of 
{X^{xi), . . .,X^{xm)) is given by 

m „ m 

E(exp(^^t,X^(x,))) =exp(-co / jj^t,- (e^<^-«> - l) [^(O-^^^-^c^e), 
j=i J^'' j=i 

for any ti, . . . ,tm G M. Thus, for any c > 0, by a change of variable 7 = c^*^ in the 
integral of the right side, since ip is an homogeneous function, we get 

m m 

E(^exp(^z^tjX^(c%)j j = E(^exp(^z^tjC-^X^(xj))), 

which proves (a). Furthermore, for any h and x G M^, we have that 

X^{x + h)-X^{x)=Re [ (e^<^'«>-l)V'(0"''~'/"W^a(f^e)- 

Hence 

m 

E(^exp(^iJ2tj {X^Xj + h)- X^(xj)))) 
= expf-co / |Vt,e'<'^'^> (e*<"^'«> - 1) "z/;(0""^"''cie 

m 

= E (^exp (^i ^ tjX^ {xj ) j j , 

proving (b). □ 

Remark 4.3. In the Gaussian case, the covariance function of the random field X<^(x) 
defined by the moving average representation (j3.1jl can be computed by an argument 
similar to Proposition 8.1.4 of [21]. Let aj = E[(X^(6'))^] for any unit vector 9, 
and define r(x) and i{x) as before so that x = T{x)^i{x). Using CoroUarv 13.21 
(a) it follows that E[(X(p(x))^] = t^xY^ af^^y and then we can use the fact that 
2X^{x)X^{y) = X^{xf + X^{yf - {X^{x) - X^{y)f to conclude that 

(4.6) E [X^ix)XM] = I H^f^'^h + riyf'^ly) ' ^(^ - vf'^l.-y)] ■ 
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In the isotropic case discussed in Remark l!0 we have r(x) = and i{x) = x/||a;||, 
and a change of variables in ()3.Hl shows that = cr^ is the same for any unit vector, 
using the fact that ip{Rx) = ip{x) for any orthogonal linear transformation R in 
this case. Then (I4.6jl reduces to the familiar autocovariance function for a fractional 
Gaussian random field. A similar argument shows that the autocovariance function 
of the random field defined by the harmonizable representation p.lj) is given by 

(4.7) E [X4x)X4y)] = \ [rixf^u^^ + r{yf"ujl^^ - r{x - 

where uj^ = E[(X^(6'))^]. For the isotropic case, where 1)4.111 reduces to the harmo- 
nizable representation HI.'M for a fractional Gaussian field, we again note that is 
constant over the unit sphere. Since a mean zero Gaussian random field is determined 
by its autocovariance function, we recover the well-known fact that the moving aver- 
age and harmonizable representations of the fractional Gaussian random field differ 
by at most a constant factor. It does not seem possible to extend this argument to 
the general case of operator scaling Gaussian random fields, since it would be difficult 
to compare ag to uJq in this case. Hence there remains an interesting open question 
under which relationship between the functions ip and ip in the Gaussian case the 
moving average representation of Theorem 13.11 and the harmonizable representation 
of Theorem 14.11 are equivalent. 

Remark 4.4. Many random fields occurring in applications have Hurst indices that 
vary with coordinate [3 IH|- Consider a random field satisfying ()1.1|) . and suppose 
that the matrix E has an eigenvector e with associated real eigenvalue A. Then it 
follows from (ILljl that the stochastic process r i— s> X{re) is self-similar with 

{X{c\e)}rm = {c''X{re)}rm for all c > 0, 

so that the Hurst index of this process is H/X. If E has a basis of eigenvectors with 
distinct real eigenvalues, then the projections of this random field onto the eigenvector 
directions yield processes with different Hurst indices in each coordinate. This also 
shows that the usual methods for estimating the Hurst index, such as rescaled range 
analysis jTHI and dispersional analysis IH), can also be applied to estimate the scaling 
indices of the operator scaling random field from data, once the proper coordinates 
are established. Estimating these coordinate directions from data is an interesting 
open question. In some practical applications, these coordinates are known from 
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the problem setup. For example, in a groundwater aquifer the coordinates of the 
hydraulic conductivity field are thought to correspond to the vertical, the direction 
of horizontal mean flow, and the horizontal direction perpendicular to the mean flow 
[Zj. In fractured rock, the scaling coordinates of the transmissivity field correspond 
to the main fracture orientations, and are usually not mutually perpendicular [22|- 
Similarly, in material science, the crack fronts determine the natural coordinates |2n| . 
We caution, however, that estimating the Hurst index in the wrong (non-eigenvalue) 
coordinates is likely to be misleading, because in those directions the field is not 
self-similar. Finally, we note that the parameters E, H in are not unique. If 
(mH) holds, then we also have {X{c^'x)} ^= {cX{x)} where E' = {1/H)E, so that 
the Hurst indices of the random field are the ratio of H and the eigenvalues of E, as 
already noted. Furthermore, the exponents of an admissible function are not unique, 
because of possible symmetries, as discussed previously in Remark I2.1fll Hence the 
Hurst index of each component is really an estimate of H/ at where < ai < ■ ■ ■ < ap 
is the real spectrum of E, and these indices, as well as the coordinate system in which 
they pertain, are the same for any choice of H and E. 

We have already seen that the OSSRFs, defined by a moving average or a harmo- 
nizable representation were stochastically continuous. In the next section we show 
that in the Gaussian case a = 2 one can get Holder regularity for the sample paths. 

5. GAUSSIAN OSSRFs 

In this section, we are interested in the smoothness of the sample paths of Gaussian 
OSSRFs given by Theorem l3.1l or Theorem 14. H respectively. Moreover we compute the 
box- and the Hausdorff-dimension of the graph of OSSRFs in these cases. We follow 
the terminology used in |H]. Using their definition of the Holder critical exponent of 
a random process (Definition 5) we state the following definition. 

Definition 5.1. Let 7 G (0,1). A random field {X{x)}x(zRd is said to have Holder 
critical exponent 7 whenever it satisfies the following two properties: 

(a) for any s G (0,7), the sample paths of X satisfy almost surely a uniform 
Holder condition of order s on any compact set, that is for any compact set 
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K gW'-, there exists a positive random variable A such that 

\X{x)-X{y)\ <A\\x-y\\' foi all x,yeK. 

(b) for any s G (7,1), almost surely the sample paths of X fail to satisfy any 
uniform Holder condition of order s. 

For a Gaussian random field X a well known result links the Holder regularity of 
the sample paths x i— > X{x,u) to those of the quadratic mean. Let us recall this 
property when the field also has stationary increments. We refer to [2] Theorem 8.3.2 
and Theorem 3.3.2 for a detailed proof. 

Proposition 5.2. Let {X{x)},j.fzRd be a Gaussian random field with stationary incre- 
ments. Let 7 G (0, 1) and assume that 

7 = sup {s > 0; E ((X(x) - X(0))') = 0||,||^o {\\xf') } ■ 

Then, for any s G (0, 7), any continuous version of X satisfies almost surely a uniform 
Holder condition of order s on any compact set. 
If moreover 

7 = inf {s > 0; ||xf ^ = oy^ii^o (E {{X{x) - X{0)f))} , 
then any continuous version of X admits 7 as the Holder critical exponent. 

The previous definition and proposition are given in [S| for random processes {d = 1) 
in order to study regularity properties of a field along straight lines. More precisely, 
when {X(x)}a.gKd is a random field, it is also interesting to study the Holder regularity 
of the process {X{xo + tu)}tm, for xq G M'' and u a unit vector. This will provide 
some additional directional regularity information. For with stationary 

increments, one only has to consider {X{tu)}te^ for all directions u. Let us recall 
Definition 6 of |H|- 

Definition 5.3. Let {X{x)}^^f^d with stationary increments and let u be any direction 
of the unit sphere. If the process {X(tu)}t£R has Holder critical exponent 7(m) we 
say that X admits 7(u) as directional regularity in direction u. 

Let us investigate these properties for the Gaussian OSSRFs given by Theorem 13. II 
or Theorem 14.11 respectively. Throughout this section we fix a real d x d matrix E 
with < ai < ■ ■ ■ < Op denoting the real parts of the eigenvalues of E. Following 
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[TH] . Section 2.1, let Vi, . . . ,Vp be the spectral decomposition of M"^ with respect to 
E. For i = 1, . . . ,p, let us denote 

W, = V,(B...(BV,, 

and Wq = {0}. Observe that E\wi has ai < ■ ■ ■ < as real parts of the eigenvalues. 
As before let q = trace(i?). 

Theorem 5.4. Let </) : M'^ — [0, oo) 6e an E -homogeneous, {P, E)- admissible func- 
tion. For < H < /3 let be the moving average Gaussian OSSRF given by The- 
orem HOI Moreover let if) : ^ [0, oo) be a continuous E^ -homogeneous function 
with ip{x) > for all x 7^ 0. For < H < ai let be the harmonizable Gaussian 
OSSRF given by Theorem \4.1\ Then any continuous version of X^ and X^,, respec- 
tively, admits H/op as Holder critical exponent. Moreover, for any i = 1, . . . ,p, for 
any directionu G Wi\Wi-i, the fields X^p and X^ admit H/ai as directional regularity 
in direction u. 

Proof. Let x G M'^. With a little abuse of notation we write X^p/^ to indicate that 
we either consider X^ or X^. Observe that X^/^{0) = and in order to apply 
Proposition 15.21 we define 

Using polar coordinates with respect to E, it is straightforward to see that 

(5.1) rj/^(x) = r(x)2^r^/4/(x)), 
where for all 9 E Sq, 

(5.2) < m < rl/^id) < M, 

since T^^^ is continuous and positive on the compact set 5*0. 

For any i = 1, ... ,p let us fix m G Wi\ Wi^i. Since the spaces Vi, . . . ,Vp are 
-E-invariant and the real parts of the eigenvalues of E\wi are ai < ■ ■ ■ < Oj it follows 
as in the proof of Lemma ED by considering the space Wi instead of M'^, that for any 
small S > there exists a constant C2 = C2{u) > such that T(tu) < C2\t\^^"''~^ for 
any \t\ < 1. Furthermore, observe that if we write u = Ui + Ui_i with n.j G Vi and 
Ui-i G Wi-i we have Ui 7^ 0. Writing tu = T{tu)^l{tu) and l{tu) = Uitu) + li-i{tu) 
with li{tu) G Vi and li-i(tu) G it follows from the i?-invariance of the spectral 
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decomposition that tUi = T{tu)^li{tu) with li{tu) ^ 0. Since we have E = Ei(B- ■ -(BEp 
where every real part of the eigenvalues of Ei equals we conclude 

= \\T{tu)'^li{tu)\\ = \\T{tu)'^H,{tu)\\ < \\T{tuf^\\\\li{tu)\\ < Cr(tM)'^'-* 

for any \t\ < 1 using the fact that ||/j(tn)|| < C3 for any \t\ < 1 and some C3 > 0. 
Hence there exists a constant Ci = Ci{u) > such that r(tn) > for any 

|t| < 1. Therefore we have shown that for all directions u & Wi \ Wi-i and any small 
6 > there exist constants 6*1,6*2 > 0, such that 

(5.3) < r(tn) < Caltr/"'"^ for all |t| < 1. 

In view of 1)5.11) . ()5.2|1 and (|5.3p we therefore get that for any direction u E Wi \ Wi-i 
and any 6 > there exist constants Ci,C2 > such that Ci|tp^/'''+'^ < r^/^(tM) < 
C2\t\'^^^°'^~^ for \t\ < 1, which by Proposition 15.21 shows that X^/^ admits H/ai as 
directional regularity in direction u. 

It follows from this that for any s G (if/op, 1) almost surely the sample paths of 
X^/^ fail to satisfy any uniform Holder condition of order s, since H/ap is the Holder 
critical exponent of X^/^, in any direction of Wp \ Wp_i. Finally, in view of 1)5.11) . 
()5.2j] and Lemma l2?T] we know that for any 6 > there exists a constant C > such 
that r^/^(x) < CWxW^^^"-""^ for < 1 and hence by Proposition [O it follows that 
any continuous version of X^/^ satisfies almost surely a uniform Holder condition of 
order s < H/ap on any compact set. This concludes the proof. □ 

Having described the Holder regularity of Gaussian OSSRFs, a natural question 
that arises is to determine the box- and the Hausdorff-dimensions of their graphs on 
a compact set. We refer to Falconer |T2] for the definitions and properties of box- 
and the Hausdorff-dimension. Let us fix a compact set K C M''. For a random field 
X on M.'^ we consider Q{X){uj) = {{x,X{x){lj))]x G K} the graph of a realization 
of this field over the compact K. We will denote dim-uGiX), resp dimsGiX), the 
Hausdorff-dimension and the box-dimension of ^(X), respectively. 

It is a well understood fact that directional regularity implies information about 
the Hausdorff-dimension of the field in that direction. See e.g. [2], Chapter 8. As an 
immediate corollary to Theorem 15.41 we get: 
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Corollary 5.5. Under the assumptions of Theorem\5.J\ we have for all i = 1, . . . ,p 



and all directions u & Wi \ Wi^i that 

dim7^{(t,X^/^(tu)) : t G [0, 1]} = 2 - H/a, a.s. 

Proof. The result is a direct consequence of Theorem 15.41 and the Corollary on page 
204 of [2j, using the fact that t t— > X^i^itu) is a /? = iY/aj-index process and that 
2 -/3 < for < /3 < 1. □ 

Our next result investigates the global box- and Hausdorff-dimension of Gaussian 
OSSRFs. 

Theorem 5.6. Under the assumptions of Theorem \5.J\ for any continuous version 
of and X^, almost surely 

dim^ Q{X^I^) = dime^(X^/^) = d+1 - H/ap. 

Proof. Let us choose a continuous version of X^/^. From Theorem 15 .41 for any s < 
H/ttp, the sample paths of X^/^ satisfy almost surely a uniform Holder condition of 
order s on K. Thus by a (i-dimensional version of Corollary 11.2 of [12], we have 



dimnG{X^/^) < dimB^(X^/^) < d+1- s, a.s. 



where dimg denotes the upper box-dimension. Therefore 



dimnG{X^/^) < dimf;^(X<^/^) < d + 1 - H/ap, a.s. 

and it remains to show that a.s. dim-nQ^X^/^) > d + 1 — H/ap. Since the lower box 
dimension satisfies diuisG{X^/^) > diuifiQ^X^/^) the proof is then complete. 
We follow the same kind of ideas developed in ^ and Let s > 1. Following the 
same argument as in Theorem 16.2 of ^2], in view of the Frostman criterion (Theorem 
4.13 (a) in [l2]), if one proves that the integral Is 

{{X^/4x) - X^/4y)f + \\x - y\\^y"^^^ 



E 

KxK 



dx dy, 



is finite, then almost surely dim>^^(X<^/^) > s. 

As before, let Vi, . . . , denote the spectral decomposition of R*^ with respect to E 
and let Wi = Vi + ■■■ + Vi. We will choose an inner-product (■, ■) on M'^ which makes 
these spaces mutually orthogonal and use the norm = (x, x)^/^. Since all norms 
on are equivalent, this entails no loss of generality. 
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Since by assumption s > 1, the function (,^^ + 1) is in L^(R) and its Fourier 



transform, denoted by fs, is not only in L°°(M) but also in L^(]R). Then we can write, 
using Fourier-inversion (fundamental lemma in (Hj) 

1 

2^ 



ie + 1] 



-s/2 



It follows that 



E 



X) 



27r" ^' 



-s/2' 



it- 



fs{t)dt 



2tt 



\x - y\ 



f2 ^{x^/4,(.^)-X^/^(.y)f) 



since X^/^ is Gaussian. Then, as fs G L"^ 



E 



fs{t)dt, 

one can find C > such that 

2 , ll„ „ l|2\-''/2' 



y\ 



E 



-1/2 



< Cm ^llx — 



y\ 



T[X 



according to and (15 .21) and using the fact that X^j^ has stationary increments. 

Let us choose A > such that K d {x E M'^; ||x|| < A/2}. Then for some constant 
C > 



I.<C 



\x\ 



r X 



-H 



dx. 



x\\<A 



as long as the integral in the right hand side is bounded. 

If p = 1, by LemmaEU for 5 > 0, one can find C > such that, for < A, 

-H ^ CWxW'^^""^'^ 



T[X 



and hence Is is finite as soon as s < d + 1 — H/ ap 



6. 



If p > 2 let us write x = Xp+y for some Xp E Vp and y G Wp^i and write x = r(x)^/(x) 
with /(x) G 5*0. Decompose l{x) = lp{x) + 9 with lp{x) G Vp and 6^ G Wp-i. By the 
direct sum decomposition we see that Xp = T{x)^lp{x) and y = t{x)^9. Moreover, 
since Vp and Wp^i are orthogonal in the chosen inner product it follows that ||a;|| < A 
implies both < A and \\y\\ < A in the associated norm. In view of the proof of 
Lemma l2T| restricted to the spaces Vp and W^p-i, respectively, it follows that for any 
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(5 > and some constants 6*1,6*2 > 0, if ||x|| < A then 

\\xp\\ < Cit{x)^^-^ and \\y\\ < C2T{xy'-^ . 
Then one can find c > such that 

t{x)^ > c\\xp\\^/''^+^ and r(x)^ > c\\y\\^/^'+^ 

and thus 

T{xf > c/2 + . 

Hence, for any 6 > 

is<c ! [ {Wx^r + \\yrY^"~'^" UKr^''^^' + Wyr^'^'^r' dyd^p- 

J\\xp\\<A J\\y\\<A 

Let k = dimVp and observe that in the present case 1 < k < d — 1. By using polar 
coordinates for both Vp and Wp„i, for some constant C* > we have Is < CJg where 

Jo Jo 

The change of variables u = tv yields 

J, = r f^^\,d~s-H/a,^S^^2 ^ -^y/2^s/2 ^^H/a,+S ^ ^,H/a,-H/a,y^ ^k^l ^^^^ 

Jo Jo 

Since — < \ < k. the second term is bounded as soon as s > k + \ — H/a„ — 5, 

dp i 

whereas the first one is finite whenever s < d + 1 — H/ap — 5. Thus, for all 5 > 
small enough, it follows that almost surely dim7^^(X^/^) > d + 1 — H/ap — 5 and 
the proof is complete. □ 

Remark 5.7. As pointed out in the introduction, the fractional Brownian sheet 
{i?//(x)}a.gigd is operator scaling with E = diag(ai, . . . , ad) where = 1/Hi, in fact 
{BH{c^x)}x(z]^d ^= {ci?//(x)}^.g]Rd, but does not have stationary increments. By The- 
orem 1.1 of [HI we know that dimT^ Q{Bh) = dimsGiBu) = (i+ 1 — min(a5~^, . . . , (y^^)- 
Now let X^/^ be our Gaussian OSSRFs considered above. It follows from operator 
scaling, that 

1 f ^ 
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too, and that = ai/H^ i = 1, . . . ,p are the real parts of the eigenvalues of {1/H)E. 
Hence Theorem 15.61 can be reformulated as 

dim-;^ g{X^/^) = dims Q{X^/^) = d + 1 - min(aj"\ . . . , ^) 

in complete similarity to the result in [2]. Hence we have constructed operator scaling 
Gaussian random fields with the same box- and Hausdorff-dimension as the fractional 
Brownian sheet, but our fields have additionally stationary increments. 
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